Astronomy & Astrophysics manuscript no. ms'0.9 


©ESO 2012 


November 21, 2012 





The young stellar population of IC 1613. III. 

New O-type stars unveiled by GTC-OSIRIS. * 

M. Garcia 1 ' 2 and A. Herrero 1 ' 2 

1 Institute de AstrofTsica de Canarias, C/ Via Lactea s/n E-38200 La Laguna, Tenerife, Spain, 
e-mail: mgg@iac.es 

2 Departamento de Astrofisica, Universidad de La Laguna, Avda. Astrofisico Francisco Sanchez, s/n, E-38071 La Laguna, Tenerife, 

<N- Spain ' 

Received july, 2012; accepted november 16th, 2012 

ABSTRACT 

>: 

Context. Very low metallicity massive stars are key to understand the reionization epoch. Radiation-driven winds are one of the 
main agents of the evolution of massive stars, and consequently an important ingredient of our models of the early-Universe. Recent 
findings hint that the winds of massive stars with poorer metallicity than the SMC may be stronger than predicted by theory. Besides 

Q\ _ calling the paradigm of radiation driven winds into question, this result would impact the calculated ionizing radiation and mechanical 

t-H . feedback of massive stars, and the role these objects play at different stages of the Universe. 

Aims. The field needs a systematic study of the winds of a large sample of very metal poor massive stars. The sampling of spectral 
types is particularly poor in the very early types. This paper's goal is to increase the list of known O-type stars in the dwarf irregular 
| galaxy IC1613, whose metallicity is smaller than the SMC's by roughly a factor 2. 

• Methods. Using the reddening-free Q-parameter, evolutionary masses and GALEX photometry, we built a list of very likely O-type 

stars. We obtained low-resolution (R~1000) GTC-OSIRIS spectra for a fraction of them and performed spectral classification, the 
only way to unequivocally confirm candidate OB-stars. 
'-M Results. We have discovered 8 new O-type stars in IC 1613, increasing the list of 7 known O-type stars in this galaxy by a factor of 

2. The best quality spectra were analyzed with the model atmosphere code FASTWIND to derive stellar parameters. We present the 
5^ ■ first spectral type - effective temperature scale for O-stars beyond the SMC. 

Conclusions. The target selection method is successful. From the pre-selected list of 13 OB star candidates, we have found 8 new 
O-stars, 4 early-B stars and provided a similar type for a formerly known early-O star. Further tests are needed but the presented 
1 1 procedure can eventually make preliminar low resolution spectroscopy to confirm candidates unnecessary. The derived effective 

temperature calibration for IC1613 is about 1000K hotter than the scale at the SMC. The analysis of an increased list of O-type stars 
t-H ' will be crucial for the studies of the winds and feedback of massive stars at all ages of the Universe. 

> " 
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in 

* . 1. Introduction Blue massive stars (BMS) experience radiation driven winds, 

powered by the scattering of photons in numerous U V transitions 
^ ; Massive stars are crucial to understand the Universe because of f metallic ions and therefore strongly dependent on metallicity. 

; their impact in many astrophysics fields. Mighty stellar winds Radiation driven winds are one of the central pillars of the cur- 
• • ■ and ionizing radiation fields, and a violent end as supernova, dis- ren t paradigm of massive stars, as the wind removes mass from 
. ^ '. ru P l their surrounding media and contaminate it with the prod- the star and changes the physical conditions at the stellar core. 

. ucts of the nuclear reactions that feed these titans. Our inter- The wind regulates, directly or indirectly the evolution of the 
^ ■ est on very low metallicity massive stars is growing rapidly be- sta r, its ionizing radiation and mechanical feedback, the super- 
03 ! cause of their role in the early Universe. In the primordial ex- nova ( SN) explosion and type of compact object left, and stellar 
tremely metal-poor environment the fo rmation of high mass yields dGeorgvet al., 2009; I Wooslev et al.L[2002b . 
stars was favoured (Broi niti et al 1 120021) and these started the 

re-ionization of the Universe. Massive stars are suspected pro- ^ ,. . . , ,. , . 

j. , , ^t,t, ur 1 — — — 'Hn Theory predicts a strong correlation between the mo- 

gemtors of long gamma ray bursts (long-GRBs Georav et al., „ J r . , , . , , , . . c 

rj rrr-, . TT n uXoA j .i • T mentum carried by the wind and the luminosity of the 

2009; Wooslev & Heaer, 2006), and this connection may ex- „ , . ... . J . , . , . ./ , .. 

' . : — , . ' — - — * tv ^r,™ -ii-i i i -j- star and metallicity: the wind-momentum luminosity relation 

plain the highest frequency of long-GRBs with higher redshirt ,„ n D i ,^ , .. , ■ / p t, , i i lr . ng |. rr, • , , .. 

F , , ? ■ , ™™ , ^ (WLR, Kudntzki, Lennon, & Puis, 1995). This relation and its 

and decreasing metallicity Gehrels et al., 2009, and references ^ n -. ~ ~, j Tr T~ ~T~ ■ . 

... ° J metallicity dependence ha ve been thouroug hly characterized 

both by theoreticians (e.g. I Vink et aTL 12001 1) and observations 

(Mokie m et all l2007t) from the Milky Way (MW) down to the 

* Based on observations made with the Gran Telescopic Canarias metallicity of the Small Magellanic Cloud (SMC). Very low 

(GTC), instaled in the Spanish Observatorio del Roque de los metallicity BMS are expected to experience weaker winds than 

Muchachos of the Instituto de Astrofisica de Canarias, in the island of SMC stars and much weaker winds than MW stars. However, 

La Palma. Program ID GTC59-1 IB. some recent results are in marked contrast. 
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We have found a resolved Luminou s Blue Variable star 
(LBV) of very low metallicity (~ 0.2Z o , iHerrero et all [2010) 
with strong optical P Cygni profiles; an analogous example ex- 
ists in NGC2366 fDris sen et all 12001) an d similar unresolved 
cases may exist in more distant galaxies dPustilnik et all 1.2008; 
llzotovetal.l 1201 ll) . rTramper et all d201 ll) reported 6 stars with 
stronger wind momentum than expected at the poor metallicity 
of their host galaxies (~ 1 /7 Z ) from X-Shooter spectroscopic 
analyses, although error bars are too large for re sults to be de- 
termi nant. Our analysis of an Of star in IC 1613 (IHerrero et all 
l20T2h concluded that the star may have a strong windor, alter- 
natively, a slower than expected wind acceleration. iLucvl d20 1 2|) 
argues that the neglect of wind clumping may explain these find- 
ings, but also suggests that the discrepancy may be caused by 
an additional wind-driving mechanism operating only at certain 
metallicities, temperat ures and luminositie s, or n egligible in the 
Galact ic counterparts. IHerrero et alj d2012l) 's and lTramper et al.l 
(2011)'s examples, if confirmed by the detailed study of a 
large sample of objects, pose a challenge to the standard the- 
ory of radiation line-driven winds, as there are few metals to 
drive the wind. Yet, this might ex plain why long-GRBs (typi- 
cally associated with type Ic SN, IWooslev & Heeer) 2006) are 
mostly found in me tal-poor environments (Mo diaz et all 12008; 
iLevesque et all 1201 Ol) but require a strong wind to remove the H 
and He envelope in the pre-SN stages. 

Besides the wind mechanical energy, the ionizing radiation 
emitted by massive stars is a chief interface of the pre-SN stages 
with the interstellar medium. The production of ionizing pho- 
tons can be estimated to first order from the effective temper- 
ature (T e ff), making T e g vs spectral type calibrations a very 
useful tool. While the temperature scale for the Milky Way 
and the Magel lanic Clouds have been addressed by a number 



of works (e.g. iGarcia & Bianchil l2004t iMartins et all |2005b; 



iMassev et alf l2009h . no characterization for lower metallicity 



environments exists. 

In order to study very metal-poor BMSs we need to reach 
out into the Loca l Group, beyond the Magellanic Clouds 
dGarcia et al. L l2oTTh. IC161 3 is the closest Local Group galaxy 
fDM= 24.27 iDolphin et all d200l ; E(B-V)=0.02 iLee et al.1 
d 1993ln with on-going star formation and poorer metallicity 
than the SMC: ~ 0.1 3Z from B-super giants (Bresolin etall 
20071). ~ 0.05Z Q from nebular studies dBresolin et all 120071 



Peimbert etafll 1 988[ iTalenl 1 1 9801: iDavidson & Kinman, 1982). 
Bresolin et al.l (120071) published a catalog of about 40 OB-stars 
in this galaxy, found from VLT-FORS2 MOS observations, but 
it only lists 6 O-types (4 of them were observed and analysed by 
Tramper et al. ( 201 ll)). These, together with the Of star analyzed 



Herrero et al. [ Mil makes a total of 7 O-type stars known in 



in 

IC1613, insufficient for a statistical characterization of the WLR 
at very low metallicities. 

It is necessary to extend the sample of known O-stars in 
IC1613, but so far low-resolution spectroscopy is needed to un- 
equivocally confirm O-candidates. In this paper we present new 
OB-stars in IC1613 found from low-resolution spectra taken 
with the Gran Telescopio Canarias (GTC). The targets were se- 
lected using a set of photometry -based criteria which works well 
for O-stars as endorsed by the derived spectral types. This work 
is part of our study of the BMS population of IC 16 13 and uses 
th e photometric cata log and study of OB ass ociations published 
in lGarcia et ail d2009l hereafter GHV09) and lGarcia et al] d20~Tol 
hereafter GHC10). All stellar identification numbers refer to the 
GHV09 catalog, unless otherwise specified. 

We present the target selection procedure optimized for O- 
type stars in Sect. [2] The observations and data reduction are 



detailed in Sect. [3] In Sect. |4]we provide spectral types for the 
targets, and in Sect.[5]we derive their stellar parameters. Finally, 
our conlusions are provided in Sect|6] 



2. Target Selection 

O-type stars cannot be chosen solely from optical photometry, 
since their colors are very similar to those of the slightly cooler 
B-types. So far spectral types are needed to unequivocally spot 
out O-stars, but this may be costly in observing time since they 
are often obscured by gas and dust and quite faint. We therefore 
need to device a method to produce a sound list of very likely 
O-type stars, and additional criteria (besides photometric colors) 
are needed to this end. 

GHV09 showed that OB-stars are found in a particular locus 
of the U-B vs Q diagram (see Fig.Q]and GHV09-Fig. 5), Q being 
the reddening-free Q pseudo-color Q = U-B - 0.72 x(B- V). 
The explanation is that Q increases monotonically towards later 
spectral types in the interval Qe[- 1.0,- 0.41, corre s pondi ng to 
O3-A0 types (see e.g. the calibration of iFitzgeraldl 0970), and 
GHC10). While Q is an indicator of spectral type, U-B holds the 
information of whether the star is reddened or not. Q is the pri- 
mary target selection criterion: starting from a catalog with small 
photometric errors (<0.05mag) to minimize the impact on Q, we 
choose "blue-Q" stars with Q < -0.8. An additional advantage 
of choosing stars mainly from their Q pseudo-color instead of 
classical B-V cuts, is that locally reddened targets are not dis- 
carded (see Fig. [2j». 

Yet, in order to separate O and B stars (see GHV09) further 
information is needed. From IC1613's "blue-Q" stars, we chose 
those with the largest evolutionary mass (> 25 M , derived by 
GHC10) that matched stellar-like sources in GALEX-FUV im- 
ages. We allowed for magnitudes as faint as V — 19.8 as long as 
the GALEX detection is clear. 

The procedure is not free of caveats. From its definition it 
follows that the Q-parameter depends on the adopted reddening 
law and is subject to larger photometric errors than individual 
filters. Besides, some O stars exhibit Q > -0.8 and would not 
be included in the current selection. Nonetheless, our selection 
criteria make the hunt for O-type stars more effective. 

The list of candidate OB-stars and their photometry is pro- 
vided in Table Q] We also included in the spectroscopic follow- 
up two objects with colors similar to the known WO and LBV- 
candidate stars in IC 1613 respective ly (see GHV09), one previ- 
ously known O star (65426, 05-6 V. lBresolin et al.l d2007l) l and 
two additional blue but fainter stars that serendipitously failed 
in-slit. Their location in IC1613 is shown in Fig. [3] Figs.Q]and[2] 
illustrate their position in the U-B vs Q and color-magnitude 
diagrams. Finding charts are provided in Appendix lAl 



3. Observations and data reduction 

The observations were carried out with the Optical System for 
Imaging and low Resolution Integrated Spectroscopy (OSIRIS) 
at the lOm-telescope GTC. The program was granted 5 hours un- 
der proposal GTC59-1 IB (P.I. M.Garcia), in service mode. The 
observations consisted on long-slit spectroscopy taken with the 
R2000B VPH and 1.2" slits. The spectra cover the -4000-5500 
A range with resolving power R ~ 1000, suitable to perform 
spectroscopic classification. The spectra are well oversampled 
(about 5 pixels) and can be rebinned to improve the signal-to- 
noise ratio (SNR). 
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Table 1. Coordinates, photometric data and association membership for the sample stars, from GHV09. Stars not belonging to any 
OB association are marked [-1] in colum ASSOC. The spectral types provided in column SpT were derived in this work. 



ID 


RA[deg] 


DEC[deg] 


ASSOC 


V 


B-V 


Q 


SpT 


Notes 




J2000.0 


J2000.0 














64066 


16.258634 


2.157802 


147 


19.03 


-0.21 


-0.97 


03 m((f» 




65426 


16.262785 


2.167927 


162 


19.62 


-0.20 


-0.91 


06III 


offslit. 05-6 V (Bresolin et al.. 2007) 


69476 


16.277063 


2.158948 


185 


18.89 


-0.22 


-0.97 


06.5 III 


36611 


16.186651 


2.109095 


44 


19.46 


-0.13 


-0.91 


07 III-V 


multiple? 


75661 


16.307566 


2.139771 


197 


19.79 


0.18 


-1.07 


08III 




67684 


16.270335 


2.159043 


175 


19.02 


-0.21 


-0.91 


08.51 


composite? 


60782 


16.249088 


2.153501 


127 


19.61 


-0.22 


-0.91 


09.5 III 


blend? 


61331 


16.250685 


2.153629 


127 


19.14 


-0.21 


-0.85 


09.7 II 


elongated shape but single detection in GHV09 


71708 


16.286754 


2.152986 


192 


19.78 


-0.02 


-1.04 


lateO III+neb 




60269 


16.247649 


2.153412 


127 


20.49 


-0.29 


-0.82 


B0.5I 


offslit 


60882 


16.249384 


2.096583 


135 


19.44 


-0.14 


-0.84 


B0.5 I-III 




27381 


16.163216 


2.146994 


-1 


18.42 


-0.07 


-0.83 


Bl-1.5 I 




35071 


16.182901 


2.084464 


-1 


18.78 


0.65 


-0.99 


B2.5 III 


V39-like colors. Double? 


44736 


16.206282 


2.106804 


56 


19.89 


-0.12 


-0.81 




offslit 


67063 


16.268118 


2.163212 


176 


19.65 


-0.04 


-0.96 






68456 


16.273197 


2.153952 


175 


21.67 


0.19 


-1.35 




WR-candidate 



Table 2. Observation log: slit position angles (PA), included targets, date of observations, seeing and transparency conditions, 
fraction of illuminated moon and airmass. 



Slit 


PA 


Stars 


Date 


Seeing 


Trans. 


Moon 


Airmass 




[deg] 






["] 




[FIT] 




slit_a 


85.6 


61331, 27381, 60269, 60782 


10-02-2011 


0.9 


Spect. 


0.40 


1.15 


slit_b 


-57.8 


69476, 71708, 75661, 65426 


10-24-2011 


1.2 


Clear 


0.05 


1.12 


slit.c 


45.85 


64066, 35071, 44736 


09-17-2011 


0.9 


Clear 


0.70 


1.24 


slit_d 


-28.6 


67684, 67063, 68456 


09-16-2011 


0.9 


Clear 


0.80 


1.44 


slit.e 


-78.8 


36611,60882 


09-17-2011 


1.0 


Clear 


0.70 


1.66" 



Morning twilight started few minutes before completing the second exposure 



0.0 



-0.5 



-1.0 - 



-1 .5 




X stars 
O WR 

+ Lote-B 
SIS Ae X A and later 

"*62024 *Tramper et i 
"*This poper 



-1 .5 



-1.0 



3.5 0.0 



0.5 



' 6 



20 



22 



x 

X* x 



w * * x 

**I &ox is* 



-0.4 -0.2 0.0 



U-E 



3.2 0.4 

B-V 



0.8 1.0 



Fig. 1. IC1613's U-B vs Q diagram. Black dots mark catalog 
stars with high quality photometry from GHV09. Other colors 
and symbols, except for filled st ars, mark the positio n of stars 
with known spectral types from iBresolin et all d2007l) . O- and 
early-B- stars concentrate in boxes 1 and 2. The pink sta r rep- 
resents 62024, an Of star anal yzed by Herrero et al. (120121) . Red 
stars represent the sample of Trampe 7et alj ( 201 1 ) in IC1613. 
Green stars mark this paper targets. 

The slits were oriented in specific angles to include several 
targets. One 1-hour long Observing Block (OB) was devoted to 
each slit. The observing conditions are compiled in Table [2] 



Fig. 2. IC1613's color-magnitude diagram. Colors and symbols 
as in Fig. [TJ OB-stars are dispersed in the color-magnitude di- 
agram, and do not concentrate on the galaxy's blue-plume. In 
sight of this diagram three targets are remarkably reddened, and 
would have been discarded with classical B-V color cuts. Note 
also that the brightest stars of the blue-plume are early-B super- 
giants. 
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Fig. 3. IC1613, GALEX's FUV-channel. North is up and East is 
left. Plus symbols mark the position of the sample stars, cho- 
sen with the following criteria: Blue Q-parameter (Q < -0.8), 
V-magnitude (V < 19.8), intense emission in GALEX images, 
large evolutionary mass (M > 25M ) and no H a emission when- 
ever possible. 



3.1. Data reduction 

Reduction was performed according to standard IRAfQ proce- 
dures. Each slit OB was split into two 1249s exposures, to enable 
cosmic ray removal. The exposures were coadded prior to any 
reduction procedure (case-A), with the imcombine routine and 
the crreject mechanism for cosmic-ray rejection. We checked 
against displacement between consecutive exposures and found 
shifts of up to 7 pixels in extreme cases (slightly larger than 
1 resolution element) in the spectral direction. For this reason 
we also reduced the individual exposures of all stars separately, 
and coadded the two spectra at the end of the reduction (case- 
B), to avoid problems with wavelength shifts. We compared the 
resulting spectra from case-A and case-B reduction and found 
undetectable differences, with a better cosmic-ray extraction in 
case-A, which we adopt. 

The ccdproc routine was used to trim the images and remove 
vignetted areas of the CCD, and also for bias subtraction and 
flat-field correction. The bias correction was performed using the 
overscan region. Each spectrum was flat-field corrected with the 
normalized attached daytime flats for that OB. After checking 
against shifts between the master arcs provided by the observa- 
tory and the attached arcs, the former were used for wavelength 
calibration because of their increased SNR. A 2-D wavelength 
calibration solution was obtained and assigned to the data using 
iraf routines identify, reidentify, fitcoords and transform. 

We note that the arc lamp image displays difuse light on the 
right part of chip-1 at A > 5500A, and the red wavelength cal- 
ibration for stars whose spectra laid in this region of the CCD 
is not reliable. The GTC observatory also reported intermitent 
technical problems with OSIRIS's active collimator during the 
dates of the observations, which may affect the stability of the 
wavelength calibration. Both factors could explain the wave- 
length calibration problem reported in Sect. [5] However, neither 
of them hamper the identification of spectral lines or the deriva- 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with the 
National Science Foundation 



tion of stellar parameters, hence have no negative impact on this 
paper goals. 

The apall task was used for target extraction and background 
subtraction. Background subtraction was complicated since a 
fraction of the targets were located in crowded regions, and some 
of them also experienced nebular contamination. The fit to back- 
ground levels was performed at the physical line of the image 
equivalent to 4500A, free from nebular transitions. If two stars 
were very close together, their apall aperture was chosen to min- 
imize contamination from the neighbor. Some targets exhibit a 
pedestal in the transversal cut, which could be caused by either 
nebular contamination or crowding. We checked H a images to 
evaluate the first option and, in the event of H„ detection, the 
sky aperture was chosen close to the star, at the pedestal. If there 
was no H a detection, then we interpreted that the plateau had 
stellar origin, and the sky contribution was chosen far from the 
pedestal. 

Finally, the R2000B VPH of GTC-OSIRIS experiences a 
ghost image, with negligible number of counts in the science or 
arc images. However, a non-negligible structure remains in the 
flat-field after normalization, which is transmitted to the science 
image after the flat-field correction. The ghost position varies 
between different flat-field images, but we have made a conser- 
vative estimate of the wavelength ranges that may be affected: 
4738-4778A (chipl) and 4736 - 4784A (chip2). 

The reduced stellar spectra are shown in Figs.|4]to|5] 



4. Spectral Classification 



Spectr al class ification was made after IWalborn & Fitzpatrickl 
jl990h 's and iLennon et al] dl992l) 's criteria for Milky Way 
stars bearin g in m ind metallicity effec ts, as explaine d in 
ICastro et all J2008) and summarized in ICastro I (1201 Oh . As 
far as we know, no spectral classification criteria exists for 
metallicity smaller than the SMC. We avoided diagnostics 
involving simultaneously metallic and non-metallic lines as 
much as possible. They may be misleading due to the dif- 
ferent metallicity of IC 1613 and the MW. Spectral sub- 
types for O-stars were derived comparing mainly the ra- 
tios of the lines HeII4541/HeI4471, HeII4200/HeI+HeII4026, 
HeII4541/HeI4387 and HeII4200/HeI4144, and luminosity 
class from the HeII4686 line and the SiIV4089/HeI4121 
ratio. The main diagnostics for B spectral subtypes were 
SiIV4089/SiIII4552, MgII4481/HeI4471, SiIII4552/SiII4128, 
SiIII4552/MgII4481 and SiII4128/HeI4121, and for luminosity 
class SiIV4089/HeI4121 and SiIII4552/HeI4387. As the lumi- 
nosity criteria often involve the undesired He-to-metal ratio, we 
checked that the observed visual magnitude, corrected from ex- 
tinction, agrees with the calibrated for the assigned luminosity 
class. The spectral classifications are provided in TableQ] 

The main sources of uncertainty are poor SNR, poorly re- 
moved or over-subtracted nebular contamination, unremoved 
cosmic -rays (exposure times were long and only 2 exposures 
per target were taken) and the difficult normalization around 
U 6 , that affects all diagnostics involving SiIV4089 and HeI4121. 
Additionally, the increased background levels of the OBs exe- 
cuted in bright time hamper the detection of the weakest spectral 
lines; it is possible that weak Hell lines are not detected in actual 
late-O stars, resulting in a bias towards B-types. 

We estimate that spectral types are uncertain by about 1-2 
sub-types, depending very much on the SNR of the spectrum and 
the nebular contamination. Nonetheless we are confident on the 
O and B classification: Hell lines are clear in the O-types, Hel 
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Fig. 4. GTC-OSIRIS spectra of OB-type stars in IC1613. The spectra have been rebinned to 3 pixels and corre cted by each star's 
radial velocity except for 71708, for which we used IC1613's systemic velocity -234 kms~ l dLu et all 1 1993b . The evolution of 
HeII4541/HeI4471 is clearly seen. 
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Fig. 5. Fig. [4] continued. 
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and Si lines are clear in the B-types, and all the spectra show stel- 
lar wings in the Balmer lines. As we have pointed out, the only 
risk in this sense is that early-B stars may actually be late-O stars 
with very weak Hell lines, undetected due to high background 
levels. 

From the maximal discovery expectation of thirteen new O- 
type stars, we have found eight O-stars and three early-B stars, 
and provided a similar spectral type for the formerly known 
early-O star 65426. 35071, whose colors where similar to the 
LBV-candidate V39, turned out to be an additional B-type star. 

The observations did not produce spectra apt for classifi- 
cation for the three remaining stars. 68456, located close to 
IC1613's WO in the U-B vs Q diagram, is dominated by neb- 
ular lines with no trace of WR's typical broad wind emissions. 
44736 and 67063 exhibit Balmer lines, and seem to also show 
Hel lines, but the spectra are too noisy for classification. 

4.1. Notes on individual targets: 

64066 (03 III((f))): The star displays strong Hell lines whereas 
Hel lines can be hardly seen. The spectrum also exhibits stellar 
Balmer wings with intense nebular lines at the core. We exam- 
ined the subtracted sky spectrum and concluded that it seems 
unlikely that the absence of Hel lines is caused by nebular con- 
tamination. Since NIII4634-40-42 is weak in emission, the star 
is classified as 03III((f)). There may be emission of NIV4058 in 
the spectrum which would render the star an ((f*)), and also of 
CIII4647-50-5 1 ; however, the SNR is too poor to be conclusive. 
64066 is very close to the 09 V star 63747, ~0.5mag fainter, but 
contamination is unlikely given the absence of Hel lines in the 
spectrum of 64066. 

65426 (06 III): This target is located inside nebulosity, and 
the [OIII] lines indicate incomplete sky subtraction. The star dis- 
plays broad stellar Balmer wings with strong nebular lines at the 
core. However, Hel lines do not seem contaminated by the neb- 
ula. HeII4686 has an artificial P Cygni-like profile caused by 
a strong cosmic-ray in one of the exposures. We only consider 
the absorption part of the line. The s tar was previously classified 
as 05-6 V bv lBresolin et al.l (2007), but we find a slightly later 
spectral type and higher luminosity class, 06 III. 

69476 (Q6 .5 III): The star is located inside GS1 
(Meaburn et all [T988I) . one of the large bubbles of the galaxy, 
and the sky subtraction is problematic. The [OIII]5007 line in- 
dicates a slight sky undersubtraction while the profile of the 
Balmer lines at the core hint that the nebular contribution has 
been oversubtracted. HeI4387 may also experience nebular con- 
tamination, but other Hel lines seem free from nebulosity. The 
Hell lines are broader than Hel. There may be emission of 
NIII4634 and CIII4647-50-51, which would render the star an 
((fc)) type, but SNR is too poor to be conclusive. 

36611 (07 III-V): This star was observed in bright time and 
almost at twilight, which explains the comparatively poor SNR 
of its spectrum. The [OIII] lines indicate severe nebular contami- 
nation, and the spectrum displays strong nebular lines at the core 
of the stellar Balmer lines. The Hell lines are broader than the 
Hel lines and seem to have substructure. The GHV09 catalog has 
two very close sources about 1" to the SW, which are fainter by 
0.8 and 2.3 magnitudes. We mark 3661 1 as possibly multiple. 

The P Cygni-like shape of HeI4471 is caused by a not to- 
tally removed cosmic ray in one of the exposures. HeII4686 
is a bit stronger than HeII4541 indicating class- V, but the 
SiIV4089/HeI4121 ratio suggests class-Ill. The star's absolute 
magnitude, if single, is also consistent with an intermediate lu- 
minosity class. 



75661 (08 III): The normalization of this star is very prob- 
lematic and hampers spectral classification. The [OIII]5007 neb- 
ular emission indicates that the nebular extraction was not com- 
plete. The spectral diagnostics disagree: the HeII4686/HeII4541 
ratio clearly indicates that the star has spectral type ear- 
lier than 08, but HeII4541/HeI4471, HeII4200/HeI4144 and 
HeII4200/HeI4026 point towards later spectral types. There 
seems to be MgII4481 in the spectrum. The weak HeII4686 
could be due to a relatively strong wind filling the line, which 
would agree with the relatively weak Hp line. The assigned type 
is a compromise of these pieces of evidence but the possibil- 
ity that the star is a binary, specially considering the difficul- 
ties in determining its radial velocity, is high. The assigned ra- 
dial velocity -250km s" 1 matches well HeII4541, HeII4686, Hg, 
HeI4922 and HeI4387, but H r , HeI4471 and HeII4200 indicate 
vrad=-350km s _1 . The luminosity class was derived extrapolat- 
ing criteria from the earlier 06 types; the star's absolute mag- 
nitude is close to the value calibrated for 08 supergiants by 
iMassevI (119981) . 

67684 (08.5 I): The [OIII] lines in absorption indicate that 
the sky was oversubtracted, which explains the strong artificial 
absorption at the core of the Balmer series. The spectral SNR is 
poor, as the star was observed at large airmass and bright time 
conditions. The normalization of the spectrum was problematic 
and hinders spectral classification; He lines are apparently found 
in broad bands with overlapped narrow lines. HeII4686 looks 
asymmetric and with several components. The spectral SNR 
prevents us from deciding whether this is an artifact from nor- 
malization or the star is multiple or actually a compact cluster. 
Nonetheless HeII4686 and HeII4200 are strong. We mark it as 
possibly composite spectrum. 

60782 (09.5 III): The star is located in a region of nebu- 
losity with intense emission in the Balmer lines. The [OIII]5007 
line indicates a slight sky oversubtraction. The spectrum may be 
contaminated by 6 1 33 1 , which is very close and 0.5 mag brigther 
but 60782's wider Balmer lines and spectral morphology in gen- 
eral suggest otherwise. There is also a nearby star with very 
similar V-magnitude and Q-color (60794, V=19.64, Q=-0.99), 
but it should be left out of the slit at the chosen position angle 
of the observations. During spectral extraction with apall addi- 
tional fainter nearby stars where detected. 

SiIV4089 is blended with H> difficulting luminosity classi- 
fication. 

61331 (09.7 II): The star is located in a region of nebulos- 
ity with intense emission in the Balmer lines. The [OIII]5007 
nebular line indicates sky oversubtraction, which explains the 
artificial absorption at the core of the Balmer lines. The appar- 
ent emission of HeII4541 is due to a cosmic-ray in one of the 
exposures. The star is very close to 60782 but contamination 
seems unlikely as the later is 0.5mag fainter. The lines of hy- 
drogen, Hel and Hell yield different radial velocities. The ap- 
parent mismatch between Hel and Hell is probably due to a 
poorly removed cosmic ray in the wing of HeI4471, and the one 
at HeII4541. The adopted radial velocity was calculated with 
Hel lines. The HeII4200 line is very weak. The luminosity indi- 
cator SiIV4089/HeI4121 suggests class between I and III, hence 
the assigned class-II. The star's absolute magnitude, calculated 
from the apparent magnitude corrected by extinction, is also con- 
sistent with an intermediate class. 

71708 (late-O III + neb.): Spectral quality is too poor to 
derive radial ve locity; IC1613's systemic velocity -234 kms~ l 
dLu et all 1 19931) is adopted. The spectrum displays broad stellar 
Balmer wings with strong nebular lines at the core. Hell lines, 
specially HeII4541 and HeII4686 are clearly seen and indicate 
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Fig. 6. Best fit FASTWIND model (red, solid line) to sample star 60782 (black). The chart shows the principal diagnostic spectral 
lines: H r and Hp (gravity), HeII4541 and HeI4471 (T eff ), and HeII4686 (log Q ). To illustrate that T eS is well constrained within 
the error bars, two additional models (non-solid lines) with T e s varying + 2000K are shown. The fit of these two models to the 
HeII4541 diagnostic line is poor, defining the error bars for temperature. 



an O-type. No spectral type can be assigned, however, since the 
spectrum displays no stellar Hel lines, very likely due to nebular 
contamination. SiIII4552 is not seen in the spectrum. HeII4686 
is strong in absorption, and stronger than HeII4541 ; according to 
the luminosity criteria, this corresponds to luminosity classes III 
and V for spectral types 08-09. Its absol ute magnitude agrees 
with the M v calibrated for a late-O giant (lMassevi.["l998l) . 

60269 (B0.5 I): The [OIII]5007 line indicates sky oversub- 
traction, which explains the asymmetry and strong absorption at 
the core of the Balmer lines. There is a velocity shift between 
the nebular lines (at IC1613's systemic velocity, measured from 
the sky spectrum) and the Balmer absorption from the final spec- 
trum. Hel and metal lines experience a smaller shift. 

This star was not a primary target of the observing run; it is 
very faint (V=20.5) and consequently its spectrum is noisy. All 
B-subtype diagnostics point towards an early-B type, but the ab- 
sence of HeI4121 hints B8-9 types. However, it is possible that 
HeI4121 is contaminated by nebular emission and this diagnos- 
tic was discarded. In absence of HeII4686 and HeI4121 lines, the 
luminosity class was assigned only from the SiIII4552/HeI4387 
ratio, resulting in class-I. However the Balmer lines are too broad 
for a supergiant: the star may be a late-O type star with un- 
detected Hell lines due to high background level compared to 
the stellar spectrum signal. Its absolute magnitude (My < -4) 
is consistent with the star being a late-O or an early-B dwarf 
star, but not with a supergiant (My = -6). The possibility 
that its spectrum is contaminated by the very nearby star 60782 
(09.5111) which is lmag brigther cannot be discarded. 

60882 (B0.5 I-III): The [OIII] lines indicate a correct back- 
ground extraction, although the profiles of the Balmer series core 
indicate a slight oversubtraction. This star was observed in bright 
time and almost at twilight, which explains the comparatively 
poor SNR of its spectrum. No Hell lines can be seen except per- 
haps HeII4541 ; the high background level of the exposures could 



be hiding the presence of weak Hell lines but the Silll triplet is 
strong and suggests a B-type. The SiIV4089/HeI4121 ratio indi- 
cates class-I but SiIII4552/HeI4387 suggest class-Ill. 

35071 (B2.5 III): The star seems double at OSIRIS ac- 
quisition image, which has enhanced resolution compared to 
GHV09's INT-WFC catalog. The criteria based on MgII4481, 
HeI4471 and SiIII4552 point towards types B2.5-3. We note, 
however, that the presence of SHI4128, expected at this spectral 
sub-type at higher metallicities, is not clear. 



5. Spectroscopic analysis 

The resolution and SNR of the dataset enabled a rough deter- 
mination of the stellar parameters of the sample O-stars. We 
used synthetic s pectra calculated with FASTWIND version 10.1 
dPuls et all l2005h that takes into account line blanketing, non- 
LTE effects and radiation driven winds, to fit mainly H r , H^, 
HeII4686, HeII4541 and HeI4471, and secondarily HeII4200, 
HeII541 1, HeI4387 and HeI4922. 

The synthetic spectra were taken from a vast grid of 
FASTWIND models with metallicit y 0.1 3Z , computed usin g 
the CONDOR facilities at the I AC dSimon-Diaz et all 1201 lh . 
0.13 Z is the equivalent met a llicity to the oxygen abun- 
dance derived by iBresolin et alj d2007l) for B-supergiants in 
IC1613. The grid covers from 28000 to 55000K, gravity ranging 
log g =2.6-4.3 (depending on temperature) and wind-strength 
Q-parameter log Q = M/(VooR*) L5 from -15.0 to -11.7. The 
exponent of the wind velocity law the helium abundance and 
the microturbulence also vary in the grid, but were kept constant 
to typical values for this analysis: /3=0.8, en e =0.09 and £=10km 
s" 1 . 

Starting from effective temperature and gravity (log g ) val- 
ues typical for the star's spectral type, we varied T e g in 1000K 
steps to fit simultaneously Hell and Hel lines. We then fixed T e g, 
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Table 3. Parameters derived for the sample stars 



ID 


SpT 


vrad 






logg 


logQ 


GHV09 




km s 


K 








64066 


03m((f)) 


-240 


> 49000 


3.8 


-12.7 


65426 


06III 


-270 


40000 


±2000 


3.8 


-12.5 


69476 


06.5111 


-350 


37000 


±2000 


3.4 


-12.5 


36611 


07III-V 


-140 


40500 


±3000 


4.3 


-15.0 


75661 


08III 


-250 


37500 


±3000 


3.9 


-12.7 


67684 


08.51 


-260 


38500 


±3000 


3.8 


-15.0 


60782 


09.5m 


-190 


33000 


±2000 


4.1 


-13.5 


61331 


09.7II 


-180 


33000 


±2000 


3.8 


-15.0 



and varied log g in 0.1 dex steps to fit the wings of the Balmer 
lines (the core was discarded since nebular contamination was 
severe, or unquantifiable, in many cases). The best-fitting grav- 
ity value usually requires the recalculation of T e g. Therefore, we 
iterated this process until all considered spectral lines except for 
HeII4686 were fitted. The observations do not include H a , but 
HeII4686 also provides information on the wind. We fitted this 
line by changing the wind strength log Q parameter. If the wind 
is not negligible it alters other lines under analysis, and the whole 
process must be iterated again to adjust temperature and gravity. 
Only when HeII4686 displays a clearly non-photospheric profile 
(e.g. asymmetric), we allowed f3 to vary. 

At the spectral resolution of the dataset, rotational and 
macroturbulent velocities could not be derived unless they are 
as high as ~250km s _1 . We also measured radial velocities for 
the targets. We noted some problems with the wavelength cal- 
ibration, specially in the red part of the spectrum, that affected 
the HeII5411 line. This cannot be explained by the shifts de- 
tected between consecutive exposures, but rather to the reported 
calibration problems in the red part of some spectra (Sect. [3]). 

The spectral quality of the dataset enabled the determina- 
tion of effective temperatures within +2000-3000 K, as shown 
in Fig. [6] Typical errors bars for log g are +0.2. We do not quote 
errors for log Q as the method provides only a rough estimate of 
the wind influence on other stellar parameters, and the derived 
values must be considered orientative. 

The spectral fit to the sample O-stars is shown in Fig. [7] 
Results are presented in table [3] 

Comments on individual targets 

64066 (03III((f))): This star displays only absorption lines 
of Hell, with no trace of Hel. We were able to derive only a 
lower limit for effective temperature. Consequently gravity is 
also poorly constrained. 

69476 (06.5111): H^ could not be used for analysis due 
to nebular contamination. The helium lines are very broad: we 
adopted vsini=250km s _1 . The HeII4541 and HeII4686 lines in- 
dicate vrad — 280km s _1 , but Hel and Balmer transitions yield 
~-350km s" 1 . 69476 is inside the bubble GS1, and the differ- 
ent radial velocities could be artificial and caused by undetected 
nebular contamination of gas masses moving at different radial 
velocity than the star. However, we examined the sky spectrum 
and there seems to be no nebular emission in Hel or Hell lines. 
Another possible explanation is that the star is a binary. We 
adopted the radial correction calculated with Hel lines. 

60782 (09.5III): The best fit to HeII4686 would actually 
require log Q =-13.0, but Balmer lines could not be fitted with 
this strong wind even if log g is increased. Contamination by 
nearby objects (see Sect. 14. II ) could be a possible explanation. 



61331 (09.7II): HeII4686 is blue-shifted and narrow, and is 
not reproduced by any combination of /3 and log Q . 

5.1. The very low metallicity temperature scale 

We present the first sub-SMC temperature scale in Fig. [8] It 
includes all O-type stars in IC1613 whose effective temper- 
ature has been derived from quantitative spectroscopic anal- 
ysis: the 8 st a rs an alyze d in this paper plus 5 stars from 
iTramper et alj (I2011I) and iHerrero et al. For compar- 

ison, we have added temp e rature det erminations o f SMC 
OB-stars from iMokiem et ail d2007l) and IMassev et all J2009), 
and the temperature scale for SMC B-supergiants derived by 
iTrundle et al.l (120071) . At the same spectral type, the temperatures 
we derive for IC1613 stars are similar to SMC stars, eventhough 
slightly hotter. 

We calculated a least squares linear regression to the temper- 
atures of the total sample of analyzed giants and supergiants in 
IC161 3. We proceeded similarly with SMC stars (Mokiem et al.l 
d2007l) and IMassev et all d2009l) ). The high end of IC1613's 
temperature scale (marked with a dotted line) was not well 
constrained, since we could only derive a lower temperature 
limit for the giant Q 3 of the sample. The other 03 star (from 
ITramper et alll201ll) is a dwarf and was not considered in the fit. 
As a test, we discarded the three stars marked as binary/multiple 
in Tableland obtained a very similar temperature scale: uncer- 
tain binarity may introduce some scatter in the derived relation 
but does not change it globally. 

The derived scale is -1000 K hotter for IC 16 13 than for the 
SMC, as expected given IC1613's poorer metal content. It may 
be argued whether the difference is significant given the error 
bars of all points considered: ~2500 K in this work vs typi- 
ca lly ~1000 K, but rea ching up to ~2000 K for some targets, 
in Mokiem et al. (2007). It is in our programmed future work to 
produce a sound temperature scale for IC 1613 from enhanced 
resolution and SNR spectra of a larger sample of stars in this 
galaxy. 



6. Summary and future work 

Very low metallicity massive stars are growing more and more 
important to our understanding of the high-redshift very metal- 
poor Universe. Current models for the first generations of stars 
rely on knowledge of present-day massive stars, with radiation 
driven winds being a main pillar. While the paradigm has been 
thouroughly tested in the Milky Way and the Magellanic Clouds, 
recent findings indicate that a different wind driving mechanism 
may be at work at poorer metallicities. The evolutionary mod- 
els for population-Ill massive stars, and the estimates of their 
ionizing, mechanical and chemical feedback would have to be 
recomputed. Only the quantitative spectroscopic analysis of a 
large sample of very metal poor massive stars can shed new light 
on these topics. 

This work uses the 10m telescope GTC as a Local Group ex- 
plorer, to unveil new OB -type stars in metal poor environments 
beyond the SMC. From low-resolution long slit spectroscopy, 8 
O-stars plus 4 early-B stars have been discovered. The total list 
of known O-type stars in IC 1613 has been increased by a factor 
of 2. A photometric -based selection method optimized towards 
O stars has been presented, which will save observing time to 
similar surveys in other galaxies. The extension of the sample 
of very metal poor O stars is a fundamental first step to test the 
theory of radiation driven winds at sub-SMC metallicity. 
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Fig. 7. Best fit FASTWIND model (red) to the sample stars (black). The chart shows the principal diagnostic spectral lines, as in 
Fig. [6] Stellar identification numbers are provided on the right-hand side of the plot. Axes are provided only for the bottom target 
for clarity sake; the same axes are used for the remaining stars. 



The GTC-OSIRIS spectra were analyzed with FASTWIND 
models to derive effective temperatures and gravities. From the 
results we have produced the first temperature scale for O-type 
stars beyond the SMC. The calibration yields higher tempera- 
tures for IC1613 stars, as expected given their comparatively 
poorer metal content, but must be used with caution given the 
large error bars of our results. 

Follow-up spectroscopy of increased resolution and en- 
hanced SNR is planned for a subset of this paper targets. Its 
quantitative analysis will provide more accurate stellar and wind 
parameters, and constraints to the wind theory. The sample will 
be as extensive as possible in its spectral type and luminosity 
class coverage so that new clues on the strong wind problems 
are found. Two interesting by-products will be produced: the first 
atlas of very metal-poor massive stars and a sound sub-SMC ef- 
fective temperature calibration. 
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Fig. A.l. Stars 60269 (B0.5 1), 60782 (09.5 III), 61331 (09.7 II), Fig. A.3. Stars 67684 (08.5 1), 69476 (06.5 III) and 71708 (late- 
64066 (03 HI((f))) and 65426 (06 III). O III+neb). 




Fig. A.2. Star 366 1 1 (07 III-V). Fig. A.4. Star 7566 1 (08 III). 



Appendix A: Finding charts for the sample stars 

Finding charts for the newly identified OB stars are provided in 
this section. The images were taken with the Wide Field Camera 
(WFC) at the 2.5m Isaac Newton Telescope (INT); for more de- 
tails see GHV09. All the identification charts show a field of 
2' x 2', observed with broad Harris V-band filter. 
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Fig. A.5. Star 60882 (B0.5 I-III). Fig. A.7. Star 3507 1 (B2.5 III). 




Fig.A.6. Star 27381 (Bl-1.5 I). 



